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We present a theory for the large suppression of the superfluid-density, p s , in BaFe 2 (Asi_ x P x ) 2 in the vicin¬ 
ity of a putative spin-density wave quantum critical point at a P-doping, x = x c . We argue that the transition 
becomes weakly first-order in the vicinity of x c , and disorder induces puddles of superconducting and antifer¬ 
romagnetic regions at short length-scales; thus the system becomes an electronic micro-emulsion. We propose 
that frustrated Josephson couplings between the superconducting grains suppress p s . In addition, the presence of 
‘normal’ quasiparticles at the interface of the frustrated Josephson junctions will give rise to a highly non-trivial 
feature in the low-frequency response in a narrow vicinity around x c . We propose a number of experiments to 
test our theory. 


Introduction. - An important focus of the study of high tem¬ 
perature superconductivity (SC) has been on the role of an¬ 
tiferromagnetism (AFM) and its relation to SC [1]. There is 
clear evidence across many different families of compounds 
that SC appears in close proximity to an AFM phase [2]; 
these families include the iron-pnictides, the electron-doped 
cuprates and the heavy-fermion superconductors. Moreover, 
the optimal transition temperature (T c ) of the SC is often sit¬ 
uated where the normal state AFM quantum critical point 
(QCP) would have been located, in the absence of supercon¬ 
ductivity. The experimental detection of the QCP is often 
challenging in the normal state, and more so in the supercon¬ 
ducting state. 

Recently, a number of measurements were reported in a 
member of the pnictide family, BaFe 2 (Asi_ x P x ) 2 , as a func¬ 
tion of the isovalent P-doping, x. The experiments show a 
phase transition involving onset of spin-density wave (SDW) 
order in the normal state above T c , which extrapolates to a 
T - 0 SDW QCP (see [3] and references therein). These ex¬ 
periments include: (/) a sharp enhancement in the effective 
mass, m*, upon approaching a critical doping from the over¬ 
doped side, as obtained from de Haas-van Alphen oscillations 
[4] and from the jump in the specific-heat at T c [5] , and, (ii) a 
vanishing Curie-Weiss temperature (i Oqw X extracted from the 
1 / T] T measurements using NMR. 

As we will review below, a number of puzzling results have 
appeared from experiments investigating whether the SDW 
QCP actually survives “under the SC dome.” Here we pro¬ 
pose a resolution of these puzzles by postulating a weakly 
first-order transition for the onset of SDW order in the pres¬ 
ence of SC order (see Fig. la). Our results are independent of 
the specific microsopic mechanism responsible for rendering 
the transition weakly first-order [6] . It is well known that ‘ran¬ 
dom bond’ disorder has a strong effect on symmetry-breaking 
first-order transitions [7], and ultimately replaces them with a 
disorder-induced second order transition in two dimensional 
systems. Our main claim is that the inhomogeneities associ¬ 
ated with these highly relevant effects of disorder can resolve 
the experimental puzzles. 

The possiblity of a QCP within the SC state was investi- 
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FIG. 1. (a) A cartoon phase-diagram showing the interplay between 
SDW and SC phases. The T N (Neel temperature) and T s (struc¬ 
tural/nematic transition) lines may survive as continuous transitions 
inside the SC phase up until the point marked V, below which they 
become weakly first-order (shown as black dashed line). The grey 
region depicts the regime where the system develops spatial inhomo¬ 
geneity due to the presence of disorder. Inset: Emulsion with SC 
grains (purple) and SDW(+SC) regions (grey), (b) The behavior of 
4 L (0) as a function of the tuning parameter, x, within different sce¬ 
narios (see text for details). 


gated by measurements [8] of the zero temperature London 
penetration depth, A L ( 0) oc 1 / y/p s (p s = superfluid-density), 
as a function of x. A sharp peak in 4^(0) was observed at 
x = x c and interpreted as evidence for a QCP [9] . However, 
this interpretation is at odds with general theoretical consid¬ 
erations [10] concerning a QCP associated with the onset of 
SDW order in the presence of a superconductor with gapped 
quasiparticle excitations [11, 12]. These considerations sug¬ 
gest that such systems will display a monotonic variation in 
4 2 (0) across the QCP, rather than a sharp peak (see dashed- 
blue/solid-red curves in Fig. lb) [10]. 

As a first step toward resolving this discrepancy, it is useful 
to place measurements of p s in the context of what is known 
about the normal state conductivity of the BaFe 2 (Asi_ x P x ) 2 
system, as these quantities are intimately related through a 
sum rule. The low temperature superfluid density of a spa¬ 
tially homogeneous superconductor can be estimated from the 
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FIG. 2. A comparison between the experimentally obtained values 
of 2^(0) [8] and those deduced from Homes’ law. The value of 2A is 
taken to be 150cm -1 [15], independent of v. 

“missing area” relation, 

2 f 2A/r 

Ps~- r I (T(z)dz, (1) 

n Jo 

where F is the elastic scattering rate and z = oj/F. In the dirty 
limit where A/F 1, the above relation yields Homes’ Law 
[13], ~ cr(0)A, whereas in the clean limit p s - p n where p n 

is the conductivity spectral weight in the normal state. Eqn. 
1 is particularly useful when the normal state resistivity data 
can reasonably be extrapolated to T = 0. By combining dc 
transport data as a function of x [14] and a measurement of 
2A from optical conductivity [15], Eq. 1 provides a lower 
bound on 2^(0) (with the assumption that A is independent of 
x). Fig.2 shows 2 2 (0) as a function of x obtained under this 
assumption (details of the procedure are presented as Supple¬ 
mentary Information). The decrease of superfluid density on 
the underdoped side reflects the growth in residual resistivity 
that begins as x drops below about 0.33. 

The values of 2 2 (0) estimated from Eq. 1 form a base¬ 
line for comparison with the experimental results presented 
in Ref. [8]. On the same graph in Fig. 2, we show the exper¬ 
imentally measured 2 2 (0) [8]. The data generally reflect the 
trend expected from the variation in the residual resistivity, 
with the exception of the sample with x = 0.3, in which the 
condensate spectral weight is suppressed by about 40% from 
the Homes’ Law estimate. Given the constraints imposed by 
the sum rule, there are two possible sources of this discrep¬ 
ancy: (/) the quasiparticle mass could be renormalized at this 
value of x, corresponding to an intrinsic decrease in p n , or, (ii) 
a considerable fraction of the (unrenormalized) p n could fail 
to contribute to the low temperature superfluid density. The 
latter possibility is suggested within the scenario that we de¬ 
velop here. 

We analyze the above experiments by assuming a weakly 
first-order transition [6], and argue that the presence of 
quenched disorder leads to formation of a micro-emulsion at 
small scales [7]. The system consists of SC puddles, where 


some of the puddles additionally have SDW order (see Fig. 
la inset). The SDW(+SC) regions, which have a locally well- 
developed antiferromagnetic moment but no long-range orien¬ 
tational order, act as barriers between the different SC grains. 
Upon moving deeper into the ordered side of the transition, the 
SDW(+SC) regions start to percolate and crossover to a state 
with long-range SDW order; this is the regime with a micro¬ 
scopically coexistent SC+SDW. As a function of decreasing 
x, the micro-emulsion is therefore a transitional state (shown 
as grey region in Fig. la) between a pure SC and a coexis¬ 
tent SC+SDW. Recent experiments in the vicinity of optimal 
doping using neutron-scattering and NMR have found results 
broadly consistent with our proposed phase diagram [16]. We 
note that the granular nature of superconductivity should have 
no effect on the bulk T c in the presence of percolating SC 
channels. 

Model. - When the system is well described in the vicinity of 
x c by a micro-emulsion as explained above, the phase fluctu¬ 
ations associated with the SC grains (shown as purple regions 
in Fig. la inset), can be modeled by the following effective 
theory, 

H 6 = ~Fj Jah COS (° a ~ 6b )> ( 2 ) 

a,b 

where J a b represent the Josephson junction (JJ) couplings be¬ 
tween grains 4 a ’ and ‘b\ We have ignored the capacitive con¬ 
tributions. 

The Josephson current across the junction will be given by 
I s = J a b sin(6 a -6b), and J a b may therefore be interpreted 
as the lattice version of the local superfluid density, p s (r ), 
i.e. J s (r) = p s {r) v s (r), with J s (r ) 9 v s (r) representing the 
superfluid-current and velocity respectively. Having a frus¬ 
trated JJ (also known as a zr—j unction) with a negative value 
of J a b leads to a local suppression in p s . Similar ideas have 
been discussed in the past in a variety of contexts (see Refs. 
[17] for a specific example), though the mechanism consid¬ 
ered here will be different. We shall now propose an explicit 
scenario under which a suppression in p s arises in the vicinity 
of putative magnetic QCPs, utilizing the SC gap structure in 
the material under question. 

The basic idea is as follows: suppose that the tunneling of 
electrons between the two grains is mediated by the SDW mo¬ 
ment in the intervening region [18], and is accompanied by a 
transfer of finite momentum that scatters them from a hole¬ 
like to an electron-like pocket. Because the SC gaps on the 
two pockets have a relative phase-difference of n, the JJ cou¬ 
pling will be frustrated [19]. 

Let us first focus on a single grain. In order to capture the 
multi-band nature of the SCs, we introduce two superconduct¬ 
ing order parameters, A* with i = ± to model the s ± state on the 
two pockets. Microscopically, these belong to regions in the 
grain having different momenta, ky, parallel to the junction. 
The gaps are related to the microscopic degrees of freedom 
[20] via the following relation, 

A ;(z) = 7 F V k||.kf l <^kf l t+k;|), (3) 

k^ 
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FIG. 3. A cartoon of a frustrated /r-junction between two supercon¬ 
ducting grains with a SDW(+SC) barrier. The SDW moment imparts 
a finite momentum transfer along the direction of the interface while 
scattering electrons from the electron (hole) pocket on one grain to 
the hole (electron) pocket on the other grain. 


where if/^ creates an electron at position z with momentum 
ky parallel to the junction and spin <x. Vk^k' is the pairing in¬ 
teraction in the Cooper channel and z is the coordinate perpen¬ 
dicular to the junction with area A. The regions % are defined 
as, R+ = {kyl&o > |k|||} and %- - {k|||&o < Ikylh where ko is an 
arbitrary momentum scale chosen such that A+ > 0, A_ < 0 
(see Fig. 3 for an illustration). We’ll assume that such a pre¬ 
scription is valid for each grain, with possibly different values 
of &o- 

Let us then write down a model for the two coupled SC 
grains with an intervening proximity coupled SDW that has a 
well developed moment, n. Our notation is as follows: we use 
a = a, b to denote the grain index and i = ± to denote the band 
index within each grain. From now on, we relabel ky as k. 
We introduce the Nambu spinor, 5ro-'^_ ko .,), 

where now creates an electron with momentum k paral¬ 
lel to the junction and at a position z (label suppressed), which 
belongs to a region of band within grain “a”. The effective 
Hamiltonian is given by, 


H g x = H a + H t , (4) 

Ha = + A -xrjv“ k(r , (5) 

a,i, k 

Ht = 8 2 n • ® ^ v 

k 

+ V a -] k A°-™' ® + H.c., (6) 


where g is the tunneling matrix element, f l (i = 0, x, y, z) act 
in Nambu space and & l (i = 0, x, y, z) act in spin space. 

In the above, H A corresponds to the bare pairing Hamilto¬ 
nian written for the ± bands within each of the two grains. H T 
represents the SDW moment mediated hopping of electrons 
from one grain to the other (represented by the a,b super¬ 
scripts) and simultaneously scattering from one band to the 
other (represented by the ± subscripts). Therefore, n imparts 
a finite momentum (along the interface) to the electrons when 
it scatters them from the electron (hole) pocket on one grain 
to the hole (electron) pocket on the other grain (shown as the 


black arrows in Fig. 3). 

Results.- Using the Ambegaokar-Baratoff relation [19], we 
can write the Josephson coupling (at T = 0) between the two 
grains as, 


Jab — 



[ Yj A e A e 

tea,Ceb 



dsz> 1 

Ei + Ei' 


(7) 


where E* = + A 2 £ and €, E represent the band indices on 

the different grains. Since AiAi' < 0, the coupling J a b < 0. 
Note that the specific nature of the frustrated tunneling arises 
from the same spin-fluctuation mediated mechanism that is 
predominantly responsible for the s ± - pairing symmetry [2]. 
However, there will also be a direct tunneling term (not in¬ 
cluded in Eqn. 4) in the Hamiltonian, which does not scatter 
the electrons from one pocket to the other, as they hop across 
the junction. The contribution to the JJ coupling from this 
term will be unfrustrated (i.e. J a b >0). 

The ratio of the tunneling amplitudes in the two different 
channels is non-universal and depends on various microscopic 
details. In particular, the emulsion is associated with a distri¬ 
bution of Josephson-couplings, P(J), with a mean coupling 
strength, (J) = J. If a substantial fraction of the JJ couplings 
become negative due to the mechanism proposed above, / will 
be small, and the superfluid density will be suppressed (see 
green curve in Fig. lb). 

We now propose a resolution as to the fate of the uncon¬ 
densed spectral weight (highlighted in Fig. 2), which can po¬ 
tentially be tested by measurements of the low frequency op¬ 
tical conductivity. Frustrated 7r-junctions host gapless states 
at the interface between the two grains [21, 22], giving rise to 
a finite density of states around zero energy (see Fig 4 inset). 
As a result of the gapless ‘normal’-fluid component at the in¬ 
terface, a fraction / of the spectral weight will be displaced 
from the superfluid-density to non-zero frequencies (shaded 
region in Fig. 4). Given that the weight of the condensate 
is proportional to /(I - /), the 40% suppression in p s for 
BaFe 2 (Asi_ JC P JC )2 in the vicinity of the putative QCP corre¬ 
sponds to / ~ 0.6. 

Our proposed optical conductivity, cr(co ), in the vicinity of 
penetration depth anomaly is shown in in Fig 4. The spectrum 
shows clearly that the connection between normal state con¬ 
ductivity and superfluid density implied by Eq. 2 will break 
down. In particular, cr dc (which is a property of the normal 
state), could vary monotonically with isovalent-doping across 
x C9 while the abundance of low-energy excitations in the im¬ 
mediate vicinity of x c would give rise to a non-monotonic 
variation in the superfluid density. This allows for an un¬ 
usual way of rearranging spectral weight in the superconduct¬ 
ing state below the gap, without violating optical sum-rules. 

The above scenario will give rise to a number of inter¬ 
esting low temperature thermodynamic and transport proper¬ 
ties, as we now discuss. First of all, there should be a strik¬ 
ing enhancement in the low-temperature thermal conductiv¬ 
ity and specific-heat, as a function of x in the narrow vicin¬ 
ity of x C9 due to the ‘normal’-component. It is important to 














4 


ct(oj) 



FIG. 4. Plot of the optical conductivity, cr(aj ), for a Drude metal 
(blue curve) and a superconductor with a large number of sub-gap 
excitations (red curve). The shaded region corresponds to the dis¬ 
placed spectral weight from the superfluid-density, p s (red arrow). 
Inset: Density of states, N(E) vs. E , for a conventional gapped su¬ 
perconductor (blue curve). A superconductor with a large number of 
subgap states has a finite N(E) below A (red curve). 


recall that this material has loop-like nodes on the electron- 
pockets [12]. However, the geometry of the electron-pockets 
and the magnitude of the gap do not change substantially in 
the vicinity of x c , and therefore it is unlikely that the contri¬ 
bution to the above quantities from the nodal-quasiparticles 
will have a drastic modificiation. It should therefore be rela¬ 
tively straightforward to disentangle the contribution arising 
from the nodal versus the ‘normal’ quasiparticles. Study¬ 
ing the NMR-spectra as a function of decreasing temperature 
(across T c ) and down to sufficiently low temperatures in the 
vicinity of x c should also reveal the spatial inhomogeneity as¬ 
sociated with the SDW regions. A large residual density of 
states in the superconducting state has been detected at a par¬ 
ticular P-doping via the power-law temperature dependence of 
1/T] ~ T [23]. Within our scenario, there should be a striking 
enhancement in this quantity as a function of doping around 
x c . Finally, we note that a promising direction for future stud¬ 
ies would be to measure the magnetic-field distribution due 
to the propagating currents in the emulsion using NV-based 
magnetometers [24]. 

Discussion. - The theoretical study in this paper was moti¬ 
vated by a number of remarkable experiments carried out in 
BaFe 2 (Asi_ x P x ) 2 , as a function of x in the normal and super¬ 
conducting phases. Our primary objective was to provide an 
explanation for the striking enhancement of the London pen¬ 
etration depth in the vicinity of a putative SDW QCP in the 
SC state. We developed a scenario based on the idea that true 
SDW criticality is masked by a weak first-order phase tran¬ 
sition in the superconducting state at T = 0. In this picture, 
quenched disorder naturally gives rise to an emulsion at small 
length scales with puddles of SC and SDW(+SC). It is then, 
in principle, possible for SDW moments at the interface of the 


SC grains to generate frustrated Josephson couplings, which 
deplete the local superfluid-density. Our proposed scenario 
naturally calls for a number of experimental tests that should 
be carried out in the near future, which should directly look for 
both the spatial inhomogeneities associated with the emulsion 
[25], and probe the gapless excitations using thermodynamic 
probes, as explained above. 

In addition to experiments on BaFe 2 (Asi_ x P x ) 2 , it should 
be important to further investigate the contrasting behavior of 
the electron-doped system, Ba(Fei_ x Co JC ) 2 As 2 , where 4^(0) 
behaves monotonically as a function of x across the puta¬ 
tive QCP [26]. Electron-doping leads to significantly higher 
^mounts of disorder compared to the isovalently-doped case, 
and would therefore lead to puddles with typically much 
smaller size [27]. Our proposed mechanism for the strong 
suppression of the superfluid-density in the isovalently-doped 
material relies on the existence of an emulsion with puddles 
of appreciable size, in the presence of an optimal amount of 
disorder. A comparison of the NMR spectra in the narrow 
vicinity of the putative QCP in the electron and isovalently 
doped materials would shed light on these microscopic differ¬ 
ences between the two families. 

Finally, though we have hypothesized that the SDW on¬ 
set transition inside the SC is, in the absence of disorder, a 
weak first order transition, we emphasize that the normal state 
properties are consistent with the presence of a “hidden” QCP 
around optimal doping [4, 5, 28]. It is plausible that in the 
normal state, different experimental techniques are probing 
the critical fluctuations associated with not one, but distinct 
QCPs as a function of x. For instance, m* extracted from 
high-field quantum oscillations is dominated by the vicinity 
of ‘hot-spots’, where quasiparticles are strongly damped due 
to coupling to the SDW fluctuations [29] . On the other hand, 
strong critical fluctuations associated with the nematic order- 
parameter [30], that couple to the entire Fermi-surface, would 
dominate m* extracted at zero-field from the jump in the spe¬ 
cific heat at T c . 
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Supplementary Material 
ESTIMATE OF A 2 L (0) FROM HOMES’ LAW 


In this supplementary material, we compare the value of the 
penetration depth obtained from experiments [8] with the pre¬ 
diction from Homes’ law; for the latter, we use a combination 
of the experimental data obtained from optical-conductivity 
and dc transport. For each value of the doping (x), we es¬ 
timate the (approximate) dc resistivity (p xx ) by extrapolating 
the curves to T = 0, from the transport data in fig. 1(b) of 
Ref. [14]. 

We estimate the value of 2A, where A is the superconduct¬ 
ing gap, from the data for optical conductivity in the super¬ 
conducting state, as shown in fig. 3(b) of Ref. [15]. Since T c 
remains relatively unchanged as a function of x in the vicinity 
of optimal doping, we assume A to be independent of x such 
that 2A « 150cm -1 (= 2.827 x 10 13 s -1 ). Then, in the dirty 
limit, 


Ps 


-<T dc A. 


( 8 ) 


In order to obtain the penetration depth, we need to restore 
various dimensionful constants such that, 



(9) 


where c(= 3 x 10 8 m/s) is the speed of light and £q(= 8.85 x 
10 -12 F/m; 1 F=1 Q _1 s) is the permitivity of free space. The 
values obtained are shown in the table below and have been 
presented in fig. 2 of the main text, along with a comparison 
to the experimental data [8]. 


X 

Pxx (from Ref. [14]) 
(ji£l cm) 

4(0) 

(pm 2 ) 

0.23 

130 

0.057 

0.27 

120 

0.053 

0.33 

35 

0.015 

0.41 

20 

0.009 

0.56 

15 

0.007 

0.64 

10 

0.004 






















